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We investigate interaction-induced valley domain walls in bilayer graphene in the zy = quantum 
Hall state, subject to a perpendicular electric field that is antisymmetric across a line in the sample. 
Such a state can be realized in a double-gated suspended sample, where the electric field changes sign 
across a line in the middle. The non-interacting energy spectrum of the ground state is characterized 
by a sharp domain wall between two valley-polarized regions. Using the Hartree-Fock approximation, 
we find that the Coulomb interaction opens a gap between the two lowest- lying states near the Fermi 
level, yielding a smooth domain wall with a kink configuration in the valley index. Our results suggest 
the possibility to visualize the domain wall via measuring the charge density difTerence between the 
two graphene layers, which we find exhibits a characteristic pattern. The width of the kink and the 
resulting pattern can be tuned by the interplay between the magnetic field and gate electric fields. 



C3 



O 
O 



> 

o 
in 

<N 

en 

o 

<N 



'i 



I. INTRODUCTION 



Two dimensional electron systems in magnetic fields 
exhibit a great richness of physics, particularly in the 
high field regime where the decreasing radius of the 
cyclotron orbits gives rise to increasing importance of 
electron-electron interactions. Two examples are the 
fractional quantum Hall effect (FQHE) and quantum 
Hall ferromagnets in the integer QHE.^'^ The essen- 
tial feature of the former is a condensation of the elec- 
trons into unusual correlated states which minimize 
the Coulomb energy, allowing the electrons to avoid 
each other as much as possible. Similarly, for the lat- 
ter. Coulomb interactions induce nonperturbative effects 
on the highly degenerate Landau bands of the non- 
interacting system. In particular, due to exchange, fer- 
romagnetism is induced in the system. A prominent 
manifestation of this state is the formation of Skyrmions 
as novel low energy excitations of the spin-polarized 
ground state (or isospin polarized states in bilayer QH 
systems). -^^ 

Quantum Hall ferromagnets have also been predicted 
for graphene in the integer quantum Hall regimes, ^^^ 
which exhibit particle- hole conjugate Landau levels and 
a peculiar i/ = QH state at zero energy. ^^^^ These two 
unique features are manifestations of the Dirac equation 
which governs the electron dynamics near the K and K' 
points in the band structure. For non-interacting elec- 
trons in graphene, four Landau levels are present near 
zero energy, associated with the two valleys and the two 
spin states. In this situation the Zeeman coupling sep- 
arates the states into two pairs above and below the 
Fermi energy. When interactions are included, the half- 
filled zero energy states spontaneously polarize due to 
exchange and give rise to a ferromagnetic ground state, ^ 
which may be spin or valley polarized depending on the 
strength of the field.— i^^ 



In addition to this interesting bulk property in the 
u = state, a coherent domain wall^i^ (DW) will be 
present between a spin polarized bulk state and an un- 
polarized region at the physical edge of a finite graphene 
ribbon. ^^'^^ This DW has also been predicted to support 
a Luttinger liquid edge mode, which is another manifes- 
tation of the Coulomb interaction in 2D systems. How- 
ever, it may be difficult to realize this spin configuration 
in currently available graphene ribbons, as their edges 
are in general rough. ^^ Moreover, such a pattern in the 
ground state is hard to probe directly. 

An "internal edge" in biased bilayer graphene (BLG) 
proposed by Martin et al. provides an alternative way 
to create a DW that circumvents the difficulty in mak- 
ing perfect physical edges. ^^'^^ This clean edge can be 
created in the middle of a bilayer graphene sample by 
placing it in an electric quadrupole gate where a poten- 
tial profile changes sign across the center of the sample, 
as shown in Fig. [TJ When the Fermi level is placed at 
zero energy, a pair of surface states with opposite chiral- 
ities and opposite isospins (valley index) are formed in 
the middle of the sample. These states are localized and 
resemble the edge states of quantum Hall systems near a 
physical edge. 

In the QH regime, bilayer graphene also exhibits 
particle-hole symmetric LLs and particle-hole degenerate 
zero energy states. Relative to the monolayer, the layer 
degrees of freedom of the bilayer system doubles the zero 
energy degeneracy. Perpendicular electric fields act as 
Zeeman fields for the layer degrees of freedom, ^^'^^ thus 
lifting their degeneracy. This effective Zeeman field can 
be tuned to be much larger than the Zeeman splitting of 
the real spin, set by the magnetic field. In the double- 
gated setting, this isospin Zeeman splitting changes sign 
in the middle of the sample, yielding level crossings simi- 
lar to the physical edge of a monolayer graphene sample. 
When interactions are included, QH ferromagnetism sets 
in and the fully polarized ground state acquires a finite 
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Figure 1: (Color online) Bilayer graphene in a double gated 
system. Upper panel: the red dot (5) represents the sublat- 
tice B on the top layer, and the black dot (A) represents the 
sublattice A on the bottom layer of a BLG. The imaginary 
line connecting the AB dimmer represents 71 bonding. The 
polarities are different on the two sides of the system. Lower 
panel: Vt marked by red solid line represents the potential 
profile for top layer, and Vb marked by black solid line rep- 
resents the potential profile for the bottom layer. Here we 
assume an adiabatic linear profile across the middle of the 
sample. 



spin stiffness. As a result, a coherent DW analogous to 
the spin DW found in Ref. 9 may form. 

In this paper, we study the interaction-induced valley 
DW in bilayer graphene in the i^ = quantum Hall state, 
in a physical configuration as shown in Fig. [TJ The per- 
pendicular magnetic field B^^ the strength of the perpen- 
dicular bias E^ , and the separation of two electric gates 
are controllable parameters. We use the Hartree-Fock 
approximation to derive the ground state and to evalu- 
ate the width of the DW in terms of these parameters. 
We find that the DW has an interlayer charge density dif- 
ference pattern, which may be accessible experimentally. 

This paper is organized as follows. In Section |TT1 we 
review the non-interacting energy spectrum of bilayer 
graphene with Bernal stacking under a perpendicular 
magnetic field and a double gated bias with different po- 
larities. In Section [ITU we derive the ground state wave- 
function and energy of the valley-kink domain wall within 
a self-consistent Hartree-Fock approximation, and eval- 
uate the size of the coherent domain wall. Section IIVI 
discusses the resulting charge density pattern. Finally, 
we summarize our results and discuss future directions 
in Section IVl 



II. NON-INTERACTING MODEL 

We consider a bilayer graphene sheet subject to a per- 
pendicular electric field which varies along the x direction 
as shown in Fig. [H With a gauge choice of A = yxB^ 
for the magnetic vector potential, the electron wavefunc- 
tions are localized in the x direction and extended along 
the y direction with a good quantum number ky. When 
the Zeeman splitting is small, the low-energy Hamilto- 
nian of biased bilayer graphene with Bernal stacking in 
the vicinity of the K valley ig^^^^'^ 
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for the Hamiltonian is 
and A, B (A, B) represent 

the sublattice wavefunctions on the top and bottom 
layer respectively. Here a = [dx -\- {x — X)] j \f2 and 
at = \—dx ^ [x — X)] />/2, where x (and all length 
scales henceforth) is in units of the magnetic length 
Ib^ the guiding center is defined as X = Is^y^ and 
ijJc = ^/2hvF/lB {vf ^ 10^ ms-^).^^^^^ The dominant 
interlayer coupling constant (71 ~ 0.3 eV) included in 
the model is between the AB dimer sites. V{x) is the 
interlayer bias, assumed to be adiabatically varied so 
that in the effective Hamiltonian for a given x, V{x) 
may be replaced by V{X). For simplicity, we consider 



V{X) = { 



X < -w/2 
{^)X, -w/2<X<w/2, 
V, X> w/2 



-V, 



(2) 



where w defines the separation of the two electric gates 
with opposite polarities, assumed to be much larger than 
Ib- The two lowest lying eigenvalues are 



and 



„K 



F^ - 



7? 



V{X) 



-V{X), 



and their corresponding eigenstates are 



(3) 



(4) 



5l,KX 



(r) = 













(5) 



p2,Kx{r) = 



AK+X)y 



( 

^o{x-X) 

'-^Mx-X) 



\ 



\ 



(6) 



X) I 



s\(X) 



s\(X) 




Figure 2: Quantum Hall energy spectrum of a double gated 
BLG, as a function of guiding center. We assume a smooth 
change in bias across x = 0, and it yields the existence of zero 
energy states at V = 0. 



where ^n{x — X) are the harmonic oscillator wavefunc- 
tions and N^ is a normalization factor. The (pi state is 
purely from the lowest Landau level (LLL) of the top 
layer, while the (j)2 state consists of the LLL from the 
bottom layer and the first LL from the top layer. Their 
distinction will be clear when we calculate the exchange 
energies. The representation of the Hamiltonian in K' 
has the same form with V{X) -^ — V(X), and a ba- 
sis where the order of components in the 4-spinor is in- 
verted, hence sf{X) = -sf{X) (with i = 1,2)^^ Note 
this means that the low-lying states of the K valley re- 
side primarily in one layer, while those of the K' valley 
are primarily in the other. 

As shown in Eqs. (|3j) and (|4]), the energy eigenval- 
ues sf^{X) and sf^ (X) are determined by the profile of 
V{X). They all vanish at X = where there is a level 
crossing between the K and K' states. We thus find 
that the noninteracting ground state of undoped BLG 
possesses two pairs of "helical" edge states with oppo- 
site chiralities from different valleys. The ground state 
is characterized by a sharp valley domain wall around 
X = 0, i.e. at zero bias. 



III. INTERACTION-INDUCED VALLEY KINK: 
HARTREE-FOCK TREATMENT 

When the Coulomb interaction is incorporated, the 
system develops a ferromagnetic nature.^ A sharp do- 
main wall between two spin states or two valleys, as 
obtained in the non-interacting ground-state described 
above, is not energetically favorable due to its large cost 
in exchange energy. The competition between the single- 
particle energy and the exchange energy gives rise to a 
lower energy state: a smooth kink in the spin/ valley de- 
grees of freedom. 



We focus on the situation shown in the center of Fig. [21 
where below the Fermi surface the filled energy states on 
the left are dominated by the K valley, and on the right 
by the K' valley. The Coulomb interaction modifies the 
sharp domain wall to a smooth valley kink which can be 
described by a trial wavefunction of the form?^ 



i*)=n 
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(7) 
Here Cl^^ and Cl^,^ create electrons in the levels e^{X) 
and e^ (X) respectively, which are closest to the Fermi 
energy, and we assume that the ei states are pushed far- 
ther away from the Fermi level when the Coulomb inter- 
action is included. |0) denotes the vacuum state where 
all lower states of negative single-particle energy are oc- 
cupied, and cp is di constant parameter. The function 
0{X) defines the valley profile of the domain wall vary- 
ing from to tt: as shown on the far left of Fig. [2l the 
filled state below the Fermi level is the K valley state 
which corresponds to ^ = 0; likewise on the far right of 
the figure, the filled state below the Fermi level is the K' 
which corresponds to ^ = tt. Eq. ([7j) may be regarded as 
a restricted Hartree-Fock approximation to the ground- 
state. In the following, we use this to study the guiding 
center dependence of 6>(X), and determine the width of 
the domain wall. 

The total Hamiltonian of the interacting electron sys- 
tem is 



H — Hq -\- Hint 5 
with single particle energy 

^0 = / ^^rxC^^CrX 



(8) 
(9) 



tX 



and interaction 



Hint = \j drdr' : p{v)V{v - r')p(r') : . (10) 

Here : (• • •) : indicates normal ordering. The density 
operator is projected into the two states closest to the 
Fermi energy, 

p{v)= J2 K'X'{r)<Prx{r)Ci,^,CrX, (H) 



,X,X' 



with (/) as defined in Eq. (|6j), r^ ^ represents the valley in- 
dex K and K^ and V(r - r^ 



—r-T — TT is the Coulomb 
interaction among the electrons {k, ~ 5.2 for a suspended 
bilayer graphene).^^'^^ We apply the Hartree-Fock ap- 
proximation to the total Hamiltonian in Eq. (|8]) and 
evaluate expectation values in the ground state given by 
Eq. ([7j). The total Hamiltonian can therefore be written 
SiS H = ^x-^x^ •> where H^^ is an effective 2x2 Hamil- 
tonian for each guiding center coordinate in the basis of 
|KX) and \K'X), 
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Jk' 



K'X.K'X 



(12) 

Here £k'x = — ^kx denote the single particle energies 
given by Eq. (|4|), and the interaction terms are 



4,,,CK<.A=sin^^, 



-^KX,KX = Eh ~ o 7 ^ (/^KX^^KX^ ) Vx,x^ (13) 
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-^K'X,K'X = £^if ~ o 7 . (^K^X^^K^X^ ) Vx,X^(14) 
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in which 



D*«^ 
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Cl^.C^x^ = cos^ ^^""^^ 
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cos — - — sm — - — . 



^H = Z)x' '^^X',x,x,X', where the integral Vx^.x^^x^^x^ 
is defined in Eqs. ()A1[) and ()A2IIA5[) . denotes the Hartree 
contribution to the single particle energies. Although Eh 
is formally divergent, in practice it is canceled by inter- 
actions with a uniform neutralizing background, which 
is not explicitly included in our Hamiltonian. The ex- 
change interaction matrix element Vx,x' = Vx',x,x',x is 
given by 



Vx,x' = -j—e 4 



Uo{X - X')Ko 



{X - X') 



r 



Ui{X-X')Ki 



{X - X') 



(16) 



where Vb = 2,2^1 n^ ^ ^^ ^^^ ^^^ modified Bessel func- 
tions which are localized at |X — X'| < 1 , and /7o,i de- 
note polynomial functions of (X — X') as described in 
Appendix [A] 

As shown in Eq. (p!2]) . the Coulomb interaction in- 
troduces off-diagonal exchange terms which open a gap. 



and yields a smooth domain wall as described by Eq. ([7]). 
The trial wave function |^) obeys the eigenvalue equation 

Using Eqs. (O and ([H]), this yields 
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where we define 

Ax 
A 



Jk.x,k'x, 
1 



B 



X 



Eh + - {Jkx,-kx + Jk'x,-k'x) , 
- {Jkx,kx — Jk'x,k'x) ■ 



This yields the relation 



\Ax\ = {e-Kx+Bx)tan0{X), 



and consequently 
sin6'(X) 



and 



cos6l(X) 



lA; 



Axf + {e-Kx + Bxf 



{s-KX + Bx) 



Axf + {eKX + Bxf 



(19) 



(20) 



The quantities |Ax| and Bx must be determined self- 
(18) consistently. 



Using Eq. (p^ for smO{X) in Eq. ([T5]) and replacing 
the sums over X' by integrals, we obtain two coupled gap 
equations, 
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(21) 



Bx = -^ I dX' 

An 



skx' + Bx' 



Ax'\^ + (e-KX' + Bx'f 



--Vx,x'- 



(22) 

When X is near the center (i.e. X <C i^), we can assume 
Ax' ^ Ax and Bx ^ Bx'- Since Ax (Ax ^ 0) is max- 
imum at X = 0, and Bx ^ at X = 0, an approximate 
solution to the gap equations takes the form 



where 



^x-^^l-ivX + Bx), (23) 



with /? = — 7i/wc5 and 



ri 



ill - ^g) V 
ill + w?) w ' 



which is the effective slope of the position dependent per- 
pendicular bias; 



Bx 
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(! + /(/?)) 



X, 



(24) 



with fip) = (8 + 8/32 + 3/34) / (8 + 8/32 + e/?-*) . 

Combining Eqs. (|23ll and ^M^ . we obtain 



lA; 



T]^ 



f(f3) 



X2 



(25) 



Substituting this into Eq. ()T8|) . this yields an expression 
for ^(X). In Fig. [3] we plot this for sample parameters 
as listed in the caption of the figure. The width of the 
valley domain wall may be estimated to be 

. Ao Ao 

T] V 

which is in general dependent on the ratio of the maximal 
Coulomb gap Aq (which has magnetic field dependence) 
to the applied bias and the separation of the two opposing 
polarity gates. 0{X) is almost linear in X at the center 
and curves up on the sides where the approximations 
of Ax' ~ Ax and Bx ~ Bx' are no longer valid. It 
is apparent that the width of the kink can be tuned by 
the interplay between the magnetic field and gate electric 
fields. We expect that an exact minimization of the trial 
wavefunction would yield a very similar result near the 
center of the DW, but the singularities in the slopes near 
the edges would be smoothed out. 




Figure 3: The guiding center dependence of 0. The param- 
eters are chosen as follows: V — 200 meV, Ib — 10 nm, 
It; = 20/s, 71/cJc = 3. 



IV. INTERLAYER CHARGE DENSITY 
PATTERN 

Here we propose a possible measurement to visual- 
ize the valley-kink domain wall derived in the previ- 
ous section. We start by projecting the density oper- 
ator of bilayer graphene into its four sublattices, i.e. 
X^P/x(r), where 



p{v) 



/^mW = Y^ ^lrx{^)^^^r'X'{Y)ClxCr'X 



t,t'X,X> 

in which \i represents the four sublattices ^4, 5, A, and 
B^ and ^^ represents the /ith component of ^kx defined 
in Eq. ^. 

Using Eq. ([7j) , the expectation value of the density on 
sublattice /i is 



0{X) 



E 



cos 



^KKxW+sin 



2 0{X) ^^ 



(26) 



e{x) ^^^ e{x) 



+e'^ cos ^ sin ^ [R^^,^ (r) + i?^,^^ (r)] ■ 



where 



Kr'xi^) = <^MrxW^/ir'x(r). 

The last term of Eq. (|26]) indicates interference between 
the K and K' valleys. In the valley transition region, 
an interference pattern would therefore be manifested by 
the charge density difference between the top (t) and the 
bottom (b) layer of the BLG, 

^p{x,y) = pt - ph, 

= {pA{r) + PB{r)) - iPAir) + P^(r)) ,(27) 
= po{x)^ PC Dw{x) cos (AKy^ip). (28) 



Here 



(29) 



{fS'^lix - X) - [1 - a{Xf] ^l{x - X)} , 

(30) 
with AK = K - K^ 

^_.^ 271FX 

Oi[X) = -^^— — ^7 — , 



and the normahzation factor Nh ^ \/l + /3^. The first 
term in Eq. ([28]) represents the average charge density 
difference between the top and the bottom layers, while 
the second term describes a charge density wave through 
which we see a rapid oscillation along the y direction 
with wave vector AK. The resulting interlayer charge 
density pattern is shown in Fig. 21 wherein the upper 
panel displays an intervalley interference along the y di- 
rection and a dipolar charge profile along the x direction. 
Across X = 0, the interlayer charge density pattern shows 
an interesting antisymmetric amplitude which is due to 
the switch in polarity of the potential profile. 
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V. CONCLUDING REMARKS 

We have proposed an experimental setup to realize a 
collective, smooth kink in the valley degrees of freedom 
for bilayer graphene at i^ = in the presence of a spa- 
tially antisymmetric bias field. The width of the kink 
is determined by an interplay between the magnetic field 
and the gate electric fields. We predict a potentially mea- 
surable interlayer charge density pattern to visualize this 
resulting electronic structure. According to Eq. (|3Q)) , the 
amplitude of the charge density pattern can be tuned by 
the ratio of ^ to 71. This pattern is possibly accessible 
to measurement, e.g. by an STM probe. 

The above results assume that the Zeeman splitting of 
the real spin is negligibly small compared to the maximal 
valley splitting set by the gate voltage. We note that for 
sufficiently strong magnetic fields where the real spin is 
resolved, two distinct crossing points appear in the non- 
interacting spectrum at zero energy, separated by a finite 
distance in real space. Consequently, a more complex 
double-kink pattern is expected to form in the interacting 
groundstate, which can be viewed as a pair of DW's with 
a mutual interaction which is tunable by the gate- volt age. 
This case will be studied elsewhere^ 

We conclude with speculations about the collective 
electronic transport behavior of this system. In analogy 
with what happens with a spin DW at the edge of single 
layer graphene at i^ = 0^'— we expect the DW to carry 



Figure 4: Interlayer charge density pattern in a BLG domain 
wall. Upper panel: a rapid oscillation along the y direction 
with wave vector AK = K — K^, and a dipolar profile along 
the X direction. Lower panel: at 2/ = , the interlayer charge 
density pattern pcDw{x,0) the thick blue curve is obtained 
from numerical integral. The parameters chosen here are: 
V = 200 meV, Ib = 10 nm, w = 20 Ib, 7i/cJc = 3, (/? = 0. 



valley currents which can lead to a valley QHE. Unlike 
the single layer case, the non-interacting energy spec- 
trum for the bilayer structure we consider has two pairs 
of states crossing the Fermi level, although one has much 
greater slope than the other. For long length scales, and 
for the purposes of static properties, in a first approx- 
imation one may ignore the higher energy states as we 
have done in this study. However, very close to x = 
the second pair of internal edge states will likely give 
the charge density profile further structure as they ap- 
proach zero energy. More importantly, these extra states 
crossing the Fermi energy open a second current-carrying 
channel, which will affect the transport properties of the 
system. An interesting set of questions in this regard is 
how the second channel couples to the first, in particular 
if they can be regarded as independent channels or if they 
are locked together by Coulomb interactions. Finally, we 
note that in the case where the splitting of real-spin is 



appreciable, the two coupled DWs are likely to support where up to corrections of order iV/^i)'^ 
a quasi ID collective mode characterized by a ladder-like 
dynamics. We leave these questions for future research. 
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Appendix A: Evaluation of the coulomb integrals 

Using Eqs. (J5]) and (J6j) we write the Coulomb integral 
as follows: 

v^x,.x..x3.x. = Ciy(i) + c^v^^^ + C3 (y(3) + y w) , 

(Al) with p = — 7i/a;c, and 



(A2) 



y(i) = I dvdv'e'^^'-^'^^y-y')%{x - X^)%{x' - X2)V(v - v')^q{x' - X3)^o(^ - X4), 

V(2) = /drdr'e^^^^-^^H2^-2^0<|>*(x-Xi)<l>*(x'-X2)V'(r-rO<l>i(x'-X3)<^i(x-X4), (A3) 

V^^^ = /drdr'e'(^^-^^H^-^0^*(x-Xi)^*(x'-X2)F(r-r0^i(x'-X3)^o(:^-X4), (A4) 

V^^^ = fdrdr'e'^^'-^'^(y-y)^l{x-Xi)^{x'-X2)V{r-r')^o{x'-Xs)^i{x-X^). (A5) 

I 

In particular, for Xi = X3 = X' and X2 = X4 = X example, Eq. ()A2p may be written explicitly in the form 
this yields the exchange interaction terms Vx^x'- As an 

I 



F^i) = - /dVdV'e<^-^')(^-^0 



exp < 



{-^ [{x - Xf + {x' - X'f + {x' - Xf + (x - X'f'\ } 



^J{x-x'f ^{y-y'f 



To evaluate this, we change variables to difference and 
center coordinates x = {x — x')^ y = y — y' , Xc = - — 2 — ^ 
and yc = 2 ^ ^^^ ^^^^ integrate over Xc and ?/c, to 

yii) ^ _y:i — g 5 — / ^^^y (A6 



We then use the 2D Fourier transform of the Coulomb 
potential 
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(A7) 
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{X - X') 



(A8) 



Similarly we can evaluate V^^\ V'^^', and V^"^' and 
obtain the final expression for the Coulomb integral in 
Eq. (HH), with 



Uo{X-X') = 13^ {X-X'f-i{X-X'f+8 



4/32 [(X - X'f 

Uiix - X') = p^{x- X'f h^ [(X - X'f - 2] - 4} . 



(A9) 
(AlO) 
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